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PACS 05. 10. -a - Computational methods in statistical physics and nonlinear dynamics 
PACS 89 . 20 . Hh - World Wide Web, Internet 

PACS 89 . 75 . Fb - Structures and organization in complex systems 

Abstract. - We study numerically the cascading failure problem by using artificially created 
scale-free networks and the real network structure of the power grid. The capacity for a vertex 
is assigned as a monotonically increasing function of the load (or the betweenness centrality). 
Through the use of a simple functional form with two free parameters, revealed is that it is indeed 
possible to make networks more robust while spending less cost. We suggest that our method to 
prevent cascade by protecting less vertices is particularly important for the design of more robust 
real- world networks to cascading failures. 
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The network robustness has been one of the most central 
topics in the complex network research [1]. In scale- free 
networks, the existence of hub vertices with high degrees 
has been shown to yield fragility to intentional attacks, 
while at the same time the network becomes robust to 
random failures due to the heterogeneous degree distribu- 
tion [2-5]. On the other hand, for the description of dy- 
namic processes on top of networks, it has been suggested 
that the information flow across the network is one of the 
key issues, which can be captured well by the betweenness 
centrality or the load [6]. 

Cascading failures can happen in many infrastructure 
networks, including the electrical power grid, Internet, 
road systems, and so on. At each vertex of the power 
grid, the electric power is either produced or transferred 
to other vertices, and it is possible that from some reasons 
a vertex is overloaded beyond the given capacity, which is 
the maximum electric power the vertex can handle. The 
breakdown of the heavily loaded single vertex will cause 
the redistribution of loads over the remaining vertices, 
which can trigger breakdowns of newly overloaded ver- 
tices. This process will go on until all the loads of the 
remaining vertices are below their capacities. For some 
real networks, the breakdown of a single vertex is suffi- 
cient to collapse the entire system, which is exactly what 
happened on August 14, 2003 when an initial minor distur- 
bance in Ohio triggered the largest blackout in the history 
of United States in which millions of people suffered with- 
out electricity for as long as 15 hours [7]. A number of as- 



pects of cascading failures in complex networks have been 
discussed in the literature [8-16], including the model for 
describing cascade phenomena [8] , the control and defense 
strategy against cascading failures [9, 10], the analytical 
calculation of capacity parameter [11], and the modelling 
of the real- world data [12]. In a recent paper [16], the cas- 
cade process in scale-free networks with community struc- 
ture has been investigated, and it has been found that 
a smaller modularity is easier to trigger cascade, which 
implies the importance of the modularity and community 
structure in cascading failures. 

In the research of the cascading failures, the following 
two issues are closely related to each other and of signif- 
icant interests: One is how to improve the network ro- 
bustness to cascading failures, and the other particularly 
important issue is how to design manmade networks with 
a less cost. In most circumstances, a high robustness and a 
low cost are difficult to achieve simultaneously. For exam- 
ple, while a network with more edges are more robust to 
failures, in practice, the number of edges is often limited 
by the cost to construct them. In brevity, it costs much to 
build a robust network. Very recently, Schafer et. al. pro- 
posed a new proactive measure to increase the robustness 
of heterogeneous loaded networks to cascades. By defin- 
ing the load dependent weights, the network turns to be 
more homogeneous and the total load is decreased, which 
means the investment cost is also reduced [15]. In the 
present Letter, for simplicity, we try to find a possible way 
of protecting networks based on the flow along shortest- 
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Fig. 1: The capacity c is assigned as c = X(l)l with the initial 
load I. The step function X(l) = 1 + a6(I/l max — /3) with two 
free parameters a and /3 is used in our model. For comparison, 
the curve for the Motter-Lai (ML) capacity model in Ref. [8], 
where = constant, is also shown. 



hop path, first proposed by Motter-Lai [8]. Through the 
use of our improved capacity model, we numerically exam- 
ine the cascades in scale-free networks and the electrical 
power grid network. Since for heterogeneously loaded net- 
works, overload avalanches can be triggered by the failure 
of only one of the most loaded vertices, the following re- 
sults are all based on the removal of one vertex with the 
highest load. Our results suggest that networks can indeed 
be made more robust while spending less cost. 

We first construct the Barabasi-Albert (BA) scale-free 
network [17] of the size N = 5000 with the average degree 
(k) ~ 4 to study the cascading failures. The BA network is 
characterized by the degree distribution p(k) ~ k^ 1 with 
the degree exponent 7 = 3, and it has been shown that the 
load distribution also exhibits the power-law behavior [6] , 
which means that there exist a few vertices with very large 
loads. 

The betweenness centrality for each vertex, defined as 
the total number of shortest paths passing through it, is 
used as the measure of the load and computed by using 
the efficient algorithm [18]. The capacity c v for the vertex 
v is assigned as 

c v — X(l v )l v , (1) 

where l v is the initial load without failed vertices. Al- 
though it should be possible to find, via a kind of the 
variational approach, the optimal functional form of X(l v ) 
which gives rise to the lower cost and the higher robust- 
ness (see below for the definitions of the two) we in this 
work simplify A(/„) as shown in Fig. [1] 
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where Q(x) = 0(1) for x < 0(> 0) is the Heaviside step 
function, l max — max„ l v , and we use a £ [0, oo) and [3 £ 
[0, 1] as two control parameters in the model. In Ref. [8] 
a constant A has been used (see Fig. Q] for comparison), 
which corresponds to the limiting case of (3 = with the 
identification A = 1 + a in our model. 



At the initial time t = 0, the vertex with the highest 
load is removed from the network, and then new loads 
for all other vertices are recomputed^ We then check the 
failure condition c v < l v {t) for each vertex, and remove 
all overloaded vertices to get the network at t + 1. The 
above process continues until all existing vertices fulfill the 
condition c v > l v {t), and the size of the giant component 
N' at the final stage is measured. The relative size of the 
cascading failures is conveniently captured by the ratio [8] 



N' 
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which we call the robustness from now on. For networks 
of homogeneous load distributions, the cascade does not 
happen and g ~ 1 has been observed [8]. Also for net- 
works of scale-free load distributions, one can have g « 1 
if randomly chosen vertices, instead of vertices with high 
loads, are destroyed at the initial stage [8]. 

In general, one can split, at least conceptually, the to- 
tal cost for the networks into two different types: On the 
one hand, there should be the initial construction cost to 
build a network structure, which may include e.g., the cost 
for the power transmission lines in power grids, and the 
cost proportional to the length of road in road networks. 
Another type of the cost is required to make the given 
network functioning, which can be an increasing function 
of the amount of flow and can be named as the running 
cost. For example, we need to spend more to have big- 
ger memory sizes and faster network card and so on for 
the computer server which delivers more data packets. In 
the present Letter, we assume that the network structure 
is given, (accordingly the construction cost is fixed), and 
focus only on the running cost which should be spent in 
addition to the initial construction cost. 

Without consideration of the cost to protect vertices, 
the cascading failure can be made never to happen by 
assigning extremely high values to capacities. However, 
in practice, the capacity is severely limited by cost. We 
expect the cost to protect the vertex v should be an in- 
creasing function of c v , and for convenience define the cost 
e as 
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It is to be noted that for a given value of cv., the original 
Motter-Lai (ML) capacity model in Ref. [8] has always a 
higher value of the cost than our model (see Fig. Al- 
though e = at /3 = 1, it should not be interpreted as 
a costfree situation; we have defined e only as a relative 
measure in comparison to the case of X(l) — 1 for all ver- 
tices. For a given network structure, the key quantities 
to be measured are g(a, (3) and e(a, (3), and we aim to in- 
crease g and decrease e, which will eventually provide us 



1 ln real situations of failures, the initial breakdown can happen 
at any vertex in the network. However, the eventual scale of dam- 
ages must be greater when a heavily loaded vertex is broken, and 
accordingly we in this work restrict ourselves to the worst case when 
the vertex with the highest load is initially broken. 
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Fig. 2: Cascading failures in the BA network of the size N = 
5000 and the average degree (k) w 4, triggered by the removal 
of a single vertex with the highest load. The robustness g 
and the cost e in Eqs. <j3j and ((4} are shown in (a) and (b), 
respectively, as functions of j3 at various a values [see Fig.[T]for 
a and P, the two parameters in the function A(Z) in Eq. ©]. 
(c) The relation between e and g at different a's. Compared 
with the ML model in Ref. [8], it is clearly shown that the 
network can be made more robust but with less cost. 



a way to achieve the high robustness and the low cost at 
the same time. 

In Fig.^a), we report the robustness g for the BA net- 
work of the size N = 5000 with the average degree (k) w 4 
as a function of /3 at a = 0.10, 0.15, 0.20, 0.25, 0.30, and 
1.0 (from bottom to top). As (3 increases further beyond 
the region in Fig. 0(a), the robustness g is found to de- 
crease toward zero (not shown here) , which is as expected 
since the larger /3 makes vertices with larger loads less pro- 
tected (see Fig. [1]). We also skip in Fig. [2] small values of 
j3 below approximately 0.001: If (3 < Zmin/7max> with the 
minimum load Z m i n , all vertices are given A(Z) — 1 + a, 
equivalent to the ML model corresponding to f3 = 0. It 
is shown in Fig. [^a) that for a < 0.30, g first increases 
and then decreases as (3 is increased, exhibiting a well- 
developed maximum g max at /3 = /?*. This is a partic- 
ularly interesting observation since the network becomes 



more robust (larger g) by protecting less vertices (larger 
(3). In more detail, the curve for a = 0.20 in Fig. [2a) 
shows the maximum g max « 0.62 (at (3* w 0.00133), which 
is about 3.5 times bigger than g » 0.175 (at = 0). In 
other words, the network can be made much more robust 
by assigning smaller capacities to vertices with less loads. 
For larger values of a, on the other hand, it is found that 
ffmax occurs at (3 = 0, which indicates that the above find- 
ing, i.e., possibility of making network more robust by 
protecting less vertices, does not hold, as exemplified by 
the curve for a = 1 in Fig. [2Ja). 

The above observation is closely related with Ref. [9], 
where it has been found that in order to reduce the size 
of cascades (or to have a larger g), some of less loaded 
vertices should be removed just after the initial attack. In 
reality, however, we believe that the direct application of 
this strategy of intentional breakdowns is not easy, for cas- 
cading failures usually propagate across the whole network 
very soon just after the initial breakdown. In contrast, we 
propose in this work a way to make the network better 
prepared to breakdowns, by protecting less vertices. 

In order to look at the cost benefit of protecting less 
vertices in a more careful way, we plot in Fig. [2jR) the 
cost e in Eq. (|4|) versus (3 at various values of a. As is 
expected from Fig. [TJ the cost e is shown to be a mono- 
tonically decreasing (increasing) function of (3 (a) at fixed 
a {(3). Take again the case with a = 0.20 as an exam- 
ple with e{(3*) w 0.153 and e{(3 = 0) = 0.2: It is then 
concluded that for a = 0.2 one can make the network 
3.5 (« 0.62/0.175) times more robust while spending only 
76.5% (w 0.153/0.2) of the original cost. 

In Fig. H{c), we use the same data as in Fig. [5{a) and 
(b), and show the relation between the robustness and 
the cost for a = 0.10, • • • , 0.30 from bottom to top. For 
comparison, the values (g,e) for (3 = 0, corresponding to 
the ML model, are also displayed as symbols at the end 
of curves. It is clearly shown that for a given a, one can 
achieve the higher robustness and the lower cost by tuning 
(3 toward the right-most point on each curve. We can also 
use Fig. [2fc) to choose the most efficient way to get a 
given robustness g: For example, suppose that g = 0.6 
is the required robustness. The vertical line for g = 0.6 
crosses several different curves, and one can choose the 
crossing point which has the lowest cost. 

We next study the cascading failures in the real net- 
work structure of the North American power grid of the 
size N = 4941 [19]. Although the electrical power grid 
network is a very homogeneous network in terms of the 
degree distribution, the load distribution, in a sharp con- 
trast, shows a strong heterogeneity as shown in Fig. [3l In 
other words, the degree distribution is more like an ex- 
ponential one, while the load distribution is similar to the 
power-law form. The broad load distribution can be one of 
the reasons of the fragility of the power grid to cascading 
failures [8]. 

We then apply, the same method as we used above, to 
the power grid, and obtain g and e as functions of (3 for 
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Fig. 3: The cumulative load distribution of power grid network 
P(l) in log-log scale. The inset shows the cumulative degree 
distribution P(k) of the power grid in linear-log scale. 
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Fig. 4: Cascading failures in the electrical power grid of the size 
TV = 4941. (Compare with Fig. [2] for the corresponding plots 
for the BA network.) The robustness g and the cost e versus 
j3 at various a values are shown in (a) and (b), respectively, 
while (c) is for the relation between e and g. Again, it is 
shown that one can achieve the higher robustness and the less 
cost simultaneously, by choosing the right-most point in (c). 



Fig. 5: Cascading failures in the BA network of the size N = 
5000 and the average degree (k) « 4, triggered by the removal 
of a single vertex with the highest load. Each vertex's capacity 
is disturbed with probability e for a = 0.2. The data are 
averaged over 20 runs. 



given values of a. Figure [4] for the cascading failures of 
the power grid is in parallel to Fig. [2] for the BA network: 
Fig. [4ja) for g versus /3, (b) for e versus /J, and (c) for e 
versus g. There are some quantitative differences between 
curves for the power grid and the BA network. However, 
qualitatively speaking, both networks are shown to ex- 
hibit the following common features: (i) For a given a, 
the robustness has a maximum <7 max at (3 = (3* , (ii) e is a 
monotonically decreasing function of (3 at a given a, and 
(hi) there exists a lob-like structure in the g-e plane, which 
indicates that one can make the network exhibit a higher 
robustness and a lower cost at the same time than the cor- 
responding values for the ML model. It is worth mention- 
ing that the power grid in Fig. |4] can be made to show the 
higher g and the lower e than the ML model in a broader 
region of a: Even at a = 1, the power grid can have much 
better robustness and much less cost in comparison to the 
ML model. Specifically, at a = 1.0 the ML model has 
g w 0.40 and e = 1.0 while our model can yield g « 0.73 
and e « 0.26 (at j3 « 0.00583) [see Fig. He)], which occurs 
when only 26% of vertices are given the higher capacity 
X(l) = 2, and the other remaining 74% of vertices have 
the lower capacity X(l) = 1. In other words, by assigning 
lower capacities to 74% of vertices, the network becomes 
much more robust. 

In reality, it is also interesting to observe the effect of 
noise on the dynamical process. In Ref. [20], when noise 
is introduced into the nonlinear dynamical system, it has 
been shown that noise changes the singularity at a special 
time to a statistical time distribution and shows various in- 
teresting behaviors. In the present work, we are interested 
in how the presence of noise influences the final cascading 
failure behavior within our scheme. Here, we introduce 
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effects of noise as an erroneous assignment of the capac- 
ity function. In detail, at a given error probability e, the 
vertex v is assigned the capacity c' v instead of its correct 
c v . 

c v = c v (l + r), (5) 



where r is the uniform random variable with zero mean 
(r G [—1,1]). We believe that this erroneous behavior 
is plausible in reality, since the perfect knowledge for the 
true value of the load for each vertex may not be available, 
which may cause an erroneous assignment of the capacity 
on a vertex. In the limiting case of e = 0, we recover our 
error-free results presented above. In Fig.[5l we report the 
results at a = 0.2 for the robustness g for the BA network 
as a function of (3 for different error probability e [see 
FigUJa) for comparison]. It is seen that for small e, the 
overall behavior is qualitatively the same as in Fig. [Ha), 
i.e., the existence of a well-developed robustness peak and 
gradual decrease as (3 is increased. The peak height of the 
robustness is found to decrease as e is increased, indicating 
the negative effect of the noise. An interesting observation 
in Fig. [5] is that as e becomes larger there exits a region 
of (3 in which the robustness is actually higher than the 
error-free case of e = 0. 

In summary, we have suggested a new capacity model to 
cascading failures, by improving the existing ML capacity 
model in Ref. [8] . The main idea in our model is the same 
as in existing studies: In a highly heterogeneous network 
with a broad load distribution, vertices with large loads 
should be more protected by assigning large capacities. 
Different from other studies in which the capacity is as- 
signed in proportion to the load, i.e., c = XI, we generalize 
the model so that the proportionality constant A is now 
changed to an increasing function X(l) of I. In more detail, 
we use the Heaviside step function for X(l) characterized 
by two parameters, the step height a, and the step posi- 
tion (3. By applying this capacity model to the artificial 
BA network as well as the real network of the power grid, 
we have clearly shown that it is indeed possible to make 
the network more robust, while at the same time the cost 
to assign capacities is drastically reduced. We believe that 
our suggested model to assign capacities to vertices should 
be practically useful in designing infrastructure networks 
in an economic point of view. As a final remark, it needs 
to be pointed out that the model proposed in this work 
should be considered as only the first step to find the op- 
timal functional form A(7) of the capacity as a function of 
the load. As a future work, we are planning to apply a 
sort of variational method to find the optimal functional 
form of X(l). 
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